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Abstract 

Metal ions play numerous important roles in biological systems being central to the function 
of biomolecules. In this letter we show that the absorption of X-rays by these ions leads to 
a complicated chain of ultrafast relaxation steps resulting in the complete degradation of their 
nearest environment. We conducted high quality ab initio studies on microsolvated Mg^"*" clusters 
demonstrating that ionisation of an Is-electron of Mg leads to a complicated electronic cascade 
comprising both intra- and intermolecular steps and lasting only a few hundreds femtoseconds. 
The metal cation reverts to its original charge state at the end of the cascade, while the nearest 
solvation shell becomes multiply ionised and large concentrations of radical and slow electron 
species build up in the metal’s vicinity. We conclude that such cascades involving metal ions are 
essential for understanding the radiation chemistry of solutions and radiation damage to metal 
containing biomolecules. 
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Q. 


Metals like Mg, Ca, Fe, Zn or Cu are essential for living organisms [11. They fulfil a 
number of important roles in biological systems, such as being instrumental for the catalytic 


activity of enzymes 


n n 


y, l2| or ensuring structural stability of chromosomes |3|, 1^. In this 
report we investigate what is the role of these naturally occurring metal species in radiation 
damage or, more broadly. X-ray induced photochemistry of biological systems. The X-ray 
absorption spectroscopy of transition metal complexes in solutions offers a glimpse at the 
transformations accompanying the interaction of X-rays with the metal. The metal ion may 


become oxidised, reduced or even remain unaltered 


BB; 


the particular outcome depends on 


the initial metal charge and the nature of ligands. Ample evidence from a different held, that 
of the X-ray crystallography, demonstrates that the radiation damage indicted on a sample 
is highly non-random a. . n.e.a.lo-p.otein. .Ke aite .Ke „.tal . 

one of the weak spots; it is often damaged at a dose much lower than the one at which the 
rest of the protein molecule is damaged js]- is common to consider the photoreduction 
of the metal ion as a manifestation of the damage {qI- This change in the metal’s charge is 
thought to change the active site’s geometry and correspondingly the diffraction image. 

In spite of accumulated knowledge mentioned above, details about the damage produced 
by the photoabsorption of an X-ray photon by the metal are sparse. Moreover, little is 
known about a sequence of events leading to this damage. In what follows we present a 
new mechanism based on high-quality ab initio calculations, which shows that the physico¬ 
chemical events started by the photoabsorption are concentrated on the metal species and 
its immediate vicinity causing thereby much damage to the molecules close to the metal. 
We hnd that the damage caused following the absorption of an X-ray photon by the metal 
can also be substantial in case the metal carries at the end of the process the same charge 
as before the photoabsorption. 

Absorption of an X-ray photon by a metal predominantly removes a core electron and 
initiates a cascade of relaxation processes. For typical metals the hrst step is Auger decay 
during which electrons with characteristic energies are emitted by the metal species which 
in turn accumulates a positive charge. The emitted Auger electrons propagate through 


the system causing secondary ionisations and resulting in damaging lesions 


10| . It is clear 


by now that in a biological medium further damaging relaxation processes are present. 
We expect the Auger decay to be accompanied and followed by ultrafast electronic decay 


processes involving the neighbouring molecules 


111 - 
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FIG. 1: Modes of electronic relaxation of the microsolvated cation created by core ionisation 

of (photoelectron not shown) and subsequent Auger decay. Excited states of Mg^"*" transfer 

the excess energy to the environment ionising a water neighbour (ICD process). Low lying states 
of Mg^"*" accept an electron from a water neighbour, transferring the excess energy either to the 
same water molecule (ETMD(2) process) or to a different one (ETMD(3) process). All described 
modes of interatomic electronic relaxation are followed by Coulomb explosions. 


These non-local processes come in two varieties (see FigH]) - interatomic Coulombic decay 
(ICD) driven b y energ y transfer, and electron transfer mediated decay (ETMD) driven by 


charge transfer 


12 


14l |. In ICD electronically excited metal species (M) de-excites by ionising 
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a neighbouring molecule (L) 


. Ln ^ M'?+ • Ln-I + L+ + eicD- 


In ETMD a neighbour donates an electron to the metal species, while the excess energy is 
used to ionise the donor (ETMD (2)) or another molecule (ETMD (3)) 

■ Ln —)■ • Ln-l + L‘^~^ + eETMD(2)i 

■ Ln —)■ • Ln -2 + 2L''' + eETMD{3)- 

Both ICD and ETMD lead to the ionisation of the nearest molecules and, in addition, to 
emitting a slow electron. However, in ICD the medium is only singly ionised, while the 
charge on the metal remains unaltered. In ETMD the medium is doubly ionised and the 
charge on the metal is reduced by one. The rate of energy transfer processes falls off as an 
inverse power of the distance, while for the charge transfer processes the rate behaviour is 
exponential [^, 15|. Therefore, previous work on the interatomic decay in rare-gas clusters 
showed that ICD proceeds on a timescale of tens of femtoseconds compared to the picosecond 
timescale of ETMD 1^. As a result, ETMD is usually observed if the decaying state does 
not have enough excess energy to decay by ICD [l7|, ll8[ . This was specihcally demonstrated 
to be the case for multiply charged ions produced in clusters by the Auger decay process 

Q. 

To study the complex chain of physico-chemical events following the X-ray photoabsorp¬ 
tion by metal atoms we considered a microsolvated cluster Mg^’''(H20)6. It may serve as a 
model of Mg^+ in solution and previous experience showed that such cluster models work 


well 


2n| . Indeed, the Mg^’*' cation usually accommodates six water molecules in its hrst 


solvation shell. The equilibrium Mg^+-0 distance in the cluster is 2.078 A and lies close 
to the corresponding value of 2.00-2.15 A in solution 21|. Therefore, since the interatomic 


decay processes involve overwhelmingly the nearest neighbours, the decay rates and branch¬ 
ing ratios should be also similar. More generally, the physico-chemical processes in the 
Mg^+(H 20)6 cluster may serve as a paradigm applicable to other metal complexes. 

We found (see Supplementary information for the computational details) that the binding 
energy of the magnesium Is-electron (K-edge) in Mg^+(H20)6 is 1317 eV. Removing this 
electron by an X-ray photon creates a highly energetic, electronically unstable trication. 
In the following we discuss the fate of this unstable trication. It predominantly decays by 


4 










emitting an Auger electron and populating a number of Mg"‘+(H20)6 states. The states 
of Mg'‘+ populated in the Auger decay are S']) (63%), Mg'^+(2s“^2p“^[^P]) 

(21%), Mg'^+(2s“^2p“^pP]) (9%), and Mg'‘+(2s“^[^S']) (7%). All of them are unstable with 
respect to intermolecular electronic decay in the presence of water molecules. Since the excess 
energy in microsolvated Mg^+(2p“^[^Zl,^ S']) relative to the Mg'^’''(2p“^pP]) ground state of 
the cation is about 4-9 eV, these states cannot decay by ICD. However, they efficiently decay 
by both ETMD(2) and ETMD(3) pathways. The computed ETMD lifetime is 16.0 fs, which 
is much shorter than the picosecond lifetimes found in rare-gas clusters 16|. Shorter ion - 
neighbour distances lead to larger orbital overlap greatly facilitating this electron transfer 


driven decay 13]. The ETMD(2) to ETMD(3) branching ratio is 1.0/1.6, therefore, both 
H20’'' and H20^’'' are efficiently produced in the ionisation of the medium. The hnal products 
are either Mg^+( 2 p“^pP])(H 20 ^+)(H 20)5 or Mg^+( 2 p“^pP])(H 20 +) 2 (H 20 ) 4 ; the charge on 
Mg decreases by one as the result of electron transfer from water (see FigH]). The ETMD 
electrons have energies between 11 and 26 eV for ETMD(2) and between 23 and 40 eV for 
ETMD(3). The greater delocalisation of the positive charge in the hnal state of ETMD(3) 
results in faster emitted electrons. 

Unlike the states considered above, microsolvated Mg^+(2s“^2p“^[^P]) possesses enough 
excess energy for the ionisation of the water molecules and, therefore, can decay by ICD. 
Its computed ICD lifetime is extremely short, 0.7 fs, showing that this interatomic decay 
occurs even faster than the already very fast local Auger decay on the metal (1.9 fs). We 
mention that the experimentally determined ICD lifetime of Mg^+(2s“^)[^S']) in liquid water 


is 1.0-1.5 fs 


22| . indicating that such extremely fast interatomic decay in hydrated metal 


ions is the rule rather than exception, due to the short metal-water distances and a large 
number of open ICD channels. In this state ICD will also dominate ETMD by about an 
order of magnitude due to the higher efficiency of energy transfer. The energies of the 
ICD electrons lie below 7 eV. Since ICD does not change the charge on Mg it produces 
Mg^+(2p“^)(H20’'')(H20)5 states. These are the same states of Mg^+ we considered above, 
but now with H20’'' ion in its vicinity. Can they continue decaying by ETMD? 

Model calculations with the water neighbour ionised via ICD replaced by a point charge 
(for details see Supplementary information) reveal that most ETMD(3) and ETMD(2) chan¬ 
nels remain open. The ETMD lifetime grows from 16.0 to 21.8 fs since only Eve neutral 
water neighbours are now available for the electronic decay. We also expect an interesting 
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interplay between ETMD and fragmentation of the Mg‘^+(2p cluster in a 


Coulomb explosion which typically proceeds on the femtosecond timescale 


23|. 


The Mg^+ states Mg^+(2s“^2p“^pP]) and Mg^+(2s“^[^S']) which are weakly populated in 
the Auger decay similarly exhibit complicated de-excitation pathways. They decay through 
a cascade of ICD and ETMD steps ultimately producing the Mg^’''(2p“^[^P]) cation. The 
lifetimes at all steps involved lie below 25 fs. 

We saw that the electronic decay after the core ionisation of magnesium leads to the 
formation of Mg^+(2p“^pP]) species. We found that both ETMD(2) and ETMD(3) channels 
are open for Mg^’''(2p“^pP])(H20)6. The corresponding lifetime is 15.5 fs, while the emitted 
electrons have energies of 0-6 and 5-19 eV, respectively. Following this decay step the 
Mg^+(H 20 ^+)(H 20)5 or Mg^+(H 20 +) 2 (H 20)4 species are produced with the ratio of 1.0 to 
1 . 2 . 

We wish to remark at this point the Mg^+(2p“^pP]) species can also be obtained directly 
from the initial core ionised Mg ion in the presence of water ligands. Thus, a 2p-electron of 
Mg may £11 the Is-vacancy and transfer its energy not to another 2p-electron, which would 


result in the Auger process on 


a core ICD-like process 


20 


m 


he metal, but to a valence electron on water ionising it in 


25|. This decay leads to population of Mg^+(H20+)(H20)5 


states. Since the positive charge in the latter is more delocalised than in the final states 
of the Auger decay, the electrons emitted in the core ICD-like process have larger energies 
(1205-1225 eV) compared to the Auger electrons (1111-1157 eV). What process will be more 
important? The Auger lifetime of the core ionised state was found to be 1.9 fs, while its ICD 
lifetime was 57 fs. Therefore, the Auger decay takes place in 97 % of Mg'^’''(H20)6 systems 
and core ICD-like only in 3%. This branching ratio, however, can be tilted much more in 


favour of the interatomic process for other metal ions. T 
ICD-like probability was found experimentally to be 10% 


lus in the case of Ca^’*' the core 


26|. 


From the previous description of the decay processes involving different Mg cations it 
is clear that these disparate steps can be glued together in one continuous cascade. Its 
schematic representation is shown in Figj2l where all steps discussed above can be found. 
Taking into account the decreasing number of neutral water neighbours available for inter¬ 
atomic decay at each consecutive step of the cascade we find that 90 percent of the core 
ionised states would cascade through to the final state within only 220 fs. Within this time 
window a large amount of Coulombic repulsion energy will be accumulated, due to the high 
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metal charge and short distances to the ionised water neighbours. Its release is expected to 
result in a complicated fragmentation pattern, involving Coulomb explosions and molecular 
dissociation |27[. Note, that after only 10 fs already 28 % of the core ionised states will un¬ 
dergo the Auger and the interatomic Coulombic decays. This sets off a Coulomb explosion 
leading to a modihcation of molecular geometry in the vicinity of the metal ion already at 
that short timescale. 

All the individual decay steps we discussed above for the microsolvated clusters will occur 
also in solution. The presence of additional solvation shells will further increase the efficiency 
of ICD, but their overall impact will be moderate, mainly because of their larger distance to 


the metal ion 


Q,E|. 


The effect of the polarisable medium will be manifested in stabilising 


of the positively charged ions at different stages of the cascade. Therefore, the energies of 
the emitted electrons will be somewhat different from the cluster’s case. 

This cascade may be viewed from two different angles: the one of the metal cation, and the 
other one of the nearest neighbours or ligands. As we mentioned above the core ionisation 
of Mg^"*" creates a highly metastable ion which first decays mostly in the Auger process. 
Even after the local Auger processes are over, the interatomic decay keeps going on. If for a 
specihc electronic state of the metal both ICD and ETMD channels are available, ICD will 
be predominant. It will lead to the electronic relaxation of the metal ion without changing 
its charge. An ETMD step will follow reducing the charge by one. This is illustrated by the 
decay of Mg‘^+(2s“^2p“^[^F]) Erst by ICD to Mg'’^+(2p“^) states and further by ETMD to 
Mg^’''(2p“^[^P]). Even if the ion does not have excess energy and ICD is not possible there 
can still be available ETMD channels. Therefore, the cascade of the interatomic processes 
goes on and the ion is being reduced until no interatomic electronic decay is allowed. For 
our system it leads to the surprising result that Mg ends up in the same electronic state it 
was in initially. 

From the point of view of the nearest neighbours the probability of their damage by the 
~1150 eV Auger electrons is small. However, they are directly damaged in the interatomic 
processes becoming either singly or doubly ionised. In the case of water solutions the single 
ionisation will lead to proton transfer reaction producing HO- radical, while doubly ionised 
water is expected to produce atomic oxygen 2^ . The hydroxyl radical is highly reactive and 
causes oxidative damage to other molecules present in the solution (see e.g. 30j). Atomic 
oxygen is reactive as well and is admittedly a source of further damaging species such as 
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FIG. 2: Schematic description of the electronic decay processes taking place after the core ionisation 
of the cation in aqueous solution. First, the intra-atomic Auger decay (black arrows) creates 

highly charged Mg^+ cations. Then, the energetic ions lose their excess energy by ICD. The Mg^"*" 
ions with little excess energy undergo ETMD and become reduced. The cascade continues until 
all Mg'^+ and Mg^+ ions revert to the original Mg^"*" state in a sequence of ultrafast, interatomic 
ICD (blue) and ETMD (red) steps. Estimated time for the complete cycle is 220 fs. Several 
reactive products such as electrons, water cations and dications are released mainly from the hrst 
coordination shell of Mg^+ as indicated by the side arrows. 


H2O2 


3l[ |. In addition to such direct damage to the nearest neighbours and the production 


of radicals, the interatomic processes result in the emission of electrons having energies <40 
eV. Such electrons can be resonantly captured by the molecules in the near environment 
initiating efficient bond breaking reactions j^. In total, the interatomic decay processes 






















massively degrade the molecules in the immediate vicinity of the metal species through 
multiple ionisations releasing in their course both reactive electrons and radicals. If we are 
to neglect the interatomic processes, the decay of the Is vacancy on Mg^’*' would result in one 
(Auger) electron. Taking them into account and counting the damaging particles released 
in the complete cascade one would obtain on average 2.4 slow electrons in addition to a fast 
Auger electron and 4.3 radicals per each Is vacancy. 

In the decay cascade presented above the Mg ion which absorbs the X-ray photon reverts 
back to its initial electronic state within few hundred femtoseconds. At the same time 
the nearest environment of this ion is multiply ionised implying the production of a large 
number of reactive particles (radicals and slow electrons) at the metal’s location. This 
shows that no change of the metal ion’s charge on X-ray irradiation is by no means equal 
to no damage done. On the ground of our detailed hndings one may anticipate that the 
similar cascades generally hold also for other metal ions independently whether the ion 
reverts back to its initial charge or photoreduction is observed. Clearly, these cascades will 
be accompanied by extensive damage to the surrounding molecules. In addition, strategic 
positioning of metal ions in biomolecules makes such damage particularly disruptive for the 
functions these molecules fulhl. Experimental evidence is available that irradiation of DNA 
molecules complexed with Ca^+ by X-rays at energies below and above the Ca K-edge showed 
about 30% enhancement in the induction of double strand breaks when the photon energy 
was increased across the K-edge 


Similarly, it was demonstrated that iron containing 


metallo-enzymes were more efficiently deactivated when irradiated with the X-rays above the 


iron’s K-edge j^|. We hope that the interatomic electronic cascades elucidated in this report 
will be useful in understanding the X-ray induced photochemistry and radiation damage of 
metal containing biomolecules. 


I. METHODS 


All computations were done by ab initio electronic structure methods, using triple-C level 


Dunning basis sets 


35l-l38|. The only exception was the calculation of Auger decay rates 


and product populations, where quintuple-C level basis sets were required. For the basis set 
details see the Supplementary information. 

The geometry of the Mg^'*’(H20)6 cluster was obtained by symmetry constrained optimi- 
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sation, relying on the M0ller-Plesset second order perturbation theory (MP2) implemented 


in the MOLPRO 2010.1 quantum chemistry package 


39l |. The T/j symmetry of the cluster 


All 


was chosen according to Feller et ah and represents the global minimum geometry 
electronic decay processes were studied at this cluster geometry. 

To identify open decay channels, the energies of the involved tricationic and tetracationic 
electronic states were calculated by the Algebraic Diagrammatic Construction (ADC(2)) 


scheme for the one-particle and two-particle propagator, respectively 


41 


^|. The latter 


allow the calculation of single and double ionisation potentials relatively to the Mg^’''(H 20 )„ 
electronic ground state introducing a perturbational expansion scheme of the propagator 
complete up to second order. The electron integrals and molecular orbitals serving as input 
or the propagator based computations were calculated by MOLCAS 7.4 software package 


44j |. The energies of pentacationic states were obtained by multi-reference conhguration 


interaction method including single and double excitations (DIRECT-CI) ^ implemented 
in the GAMESS-UK 8.0 quantum chemistry package The reference space was built up 
starting with the Mg^’''(H20)„ ground state Hartree-Fock determinant. The hnal states of 
ETMD were constructed by selecting all conhgurations having a hole in a 2s- or 2p-orbital 
of Mg plus two additional holes in the valence orbitals localised on water molecules. The 
hnal states of ICD were constructed by selecting all conhgurations having two holes in 2s- 
2p-orbitals of Mg plus an additional hole in the valence orbitals localised on water molecules. 

The total and partial electronic decay widths were calculated by means of Fano-ADC- 


Stieltjes method Due to numerical limitations, t 


states were determined in an additive approximation 


le decay widths for the Mg^+(H 20)6 


48 


49| . For details of this approxima¬ 


tion and calculation of the partial decay widths see the Supplementary information. 
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I. IDENTIFICATION OF THE ELECTRONIC STATES INVOLVED IN THE DE¬ 


CAY CASCADE 

Core ionisation of the Mg^’''(H20)6 cluster creates a highly energetic magnesium trication 
which decays by the Auger mechanism and accumulates positive charge. A series of in¬ 
teratomic electronic decay steps leads to the reduction of magnesium to the original Mg^+ 
state. The following steps of electronic decay were considered in this work: 


Mg^+{ls-^){H20\ iCD-iike^ Mg^+{H20+)iH20)5 


Mg^+iH^O), ^ Mg^+{H20+){H20), 

Mg^+{H20)e Mg^+{H20‘^^){H20), ( 2 ) 

Mg^^{H20)e Mg^^{H20^)2{H20), 


Mg^+{H20), ^ Mg^+{H20+){H20), 

Mg^+{H20), Mg‘^^{H20^+){H20), (3) 

Mg^+iH20)e Mg^+iH20+)2iH20),. 


The involved electronic states comprise tricationic, tetracationic and pentacationic states. 
For these states the energies had to be computed in order to determine open decay channels. 
All electronic structure calculations were carried out at the equilibrium geometry of the 
Mg^+(H 20)6 cluster using the ground electronic state of Mg^+(H 20)6 as the reference state. 
The initial Mg^’''(H20)6 state energy was determined by the Mpller-Plesset second order 
perturbation theory (MP2). 

We start by considering the tricationic states, which appear as initial states in Eqs. 
<m- Their energies were obtained by the Algebraic Diagrammatic Construction (ADC(2)) 


method: a propagator approach based on perturbational expansion of the Green’s 


unction 
The 


UJ 

Q 


complete up to second order. The ADC(2)-extended scheme was used throughout 
energies of the tetracationic states were calculated by the ADC(2) method for two-particle 
propagator 


2 












The ETMD and ICD of high lying tetracationic states lead to pentacationic states (Eq. 
()2|)). For these states no propagator approach was available. Therefore, a mnlti-reference 
conhgnration interaction method inclnding single and donble excitations was nsed. The con- 
hgnration interaction method is compntationally demanding and not size extensive. These 
drawbacks limited the size of the model system to the metal cation and two water neigh- 
bonrs. The energy difference between the initial tetracationic and the hnal pentacationic 
states determines the kinetic energy of the ICD and ETMD electrons. To estimate the error 
introdnced into the calcnlation of this energy difference by retaining only two water ligands 
instead of six we analysed the decay of the tricationic states (Eq. ([3])). The corresponding 
ETMD electron energies revealed a decrease of 1-2 eV when n was increased from 2 to 6. 
For the energies of ICD electrons an opposite trend is observed: they grow by 3-4 eV as n 
increases from 2 to 6. The energies of the ETMD electrons for Mg^+(Il20)2 were fonnd to 
lie above 10 eV, thus all ETMD channels open in Mg^+(H20)2 are expected to remain open 
in Mg'^’''(H20)6, while the ETMD electron energies are expected to be shifted by few eV 
to lower values. The ICD electrons, lying below 7 eV for Mg"^+(Il20)2, are expected to be 
shifted to higher values. 


II. DETERMINATION OF THE WIDTHS OF ELECTRONIC DECAY 


In order to discuss the timescales on which the electronic decay processes occur, the corre¬ 
sponding widths were calculated. For this purpose we used the Fano-Stieltjes-ADC approach. 
Within this ab initio approach the bound and continuum parts of the decaying state are 
constructed by means of the ADC(2)-extended scheme for the one-particle and two-particle 
propagator, respectively [4, . The renormalisation of the pseudocontinuum generated by 

the ADC scheme in an basis is done by the Stieltjes imaging technique j^. The partial 
decay widths were estimated by applying projector operators to the pseudocontinuum com¬ 
ponent of the wave function. The projector operators were dehned by a channel or a group 
of channels of interest (e.g. all ETMD(2) channels). 

The limitation of the Fano-Stieltjes-ADC method is given by the size of the electronic inter¬ 
action matrices to be diagonalised for the construction of the bound and pseudocontinuum 
components of the resonance wave function. For example, the interaction matrices of the 
tetracationic decaying states are essentially spanned by the large Ship block, making a 
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full ab initio computation of decay widths unfeasible for bigger systems. However, as was 
demonstrated in previous studies on ETMD and ICD, a good approximation to the elec¬ 
tronic dec^ widths can be obtained by adding up decay widths calculated for subunits of a 
system [3, . 



FIG. 1: Division of the Mg'^'*“(H20)6 cluster into subunits required for the calculation of the decay 
widths by the additive approach. Note: the symmetry lowering may mix the states belonging to 
the different irreducible representations of the T/j group of the Mg'^'*'(H 20)6 cluster. 


The approximation can be exemplihed for the total ICD width in a cluster with N water 
neighbours: 

N 

= (4) 

i=l 

where Tf*"^ is the ICD width due to the i-th water neighbour (see FigJT]). The total ETMD(2) 
width is obtained in a similar way: 

N 

= ( 5 ) 

i=l 

For the determination of the total ETMD(3) width the system is divided into subunits 
containing the metal ion and pairs of neighbours ij: 

= ( 6 ) 

i<j 

The pairs in the case of the Mg'?+(H20)6 clusters (see Fig. [1]) can further be divided into 
trans and cis type pairs: 
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r'^" = + E r!f (<:«)• (7) 

The Mg^+(H20)6, Mg^+(H20)6 electronic states correlate with the atomic multiplets such as 
Mg^’''(2p“^[^D]). The multiplets are weakly split (<0.1 eV) in the ligand held of the water 
molecules. Moreover, we assume that they are equally populated in the electronic decay. To 
simplify the calculations while retaining their accuracy we computed the interatomic decay 
widths averaged over such multiplets: 


r = yEr™. (8) 

m 

where I denotes the multiplicity and m an individual member of the multiplet. A total 
decay width of Mg'^+(H 20)6 is calculated by the summation of the multiplet averaged decay 
widths due to the subunits: 


r = E r. = E 7 E (9) 

i i m 

We compared the results of a full ab initio calculation with the values determined by the ad¬ 
ditive approach for the Mg^’''(2p“^pP])(H20)6 states. The corresponding multiplet averaged 
ETMD(2) width is given by: 

r‘" = E5E(rS’)i = 5Er®. (lo) 

2=1 m=l m=l 

For the ETMD(3) width we obtain: 

= ^m\trans) + y (H) 

m=l m=l 

The ETMD lifetime of Mg^’''(2p“^pP])(H20)6 obtained by the additive approach (15.6 fs) 
lies close to the value obtained by the full ab initio calculation (17.3 fs). The branching ratios 
of ETMD(3) to ETMD(2) product populations are 1.50 and 1.23, respectively. Therefore, 
we conclude that the approximate method is able to reproduce the total decay width and 
the branching ratios fairly well. 

The additive approach was utilised to compute the decay rates of the Mg^+(H20)6 states. 
We assumed that all decay channels which are open in a subunit with a single neighbour 
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TABLE I: Decay lifetimes of cations produced by Auger process in the core ionised Mg^'’'(H20)6 
cluster. Minor channels, contributing < 10 percent to the product population, are not shown. 


Decaying state Lifetime [fs] ETMD(3) ETMD(2) ICD 


2p-2[ip] 

15.6 

0.65 

0.35 


2p-2p5] 

18.2 

0.57 

0.43 


2s"^2p-ipP] 

17.5 

0.51 

0.49 


2s-^2p-ipP] 

0.7 



1.0 

2s-2pS] 

0.9 



1.0 


{n = 1) or a pair of neighbours {n = 2) will remain open in the larger {n = 6) system (see the 
discussion in the previous section). In the Table Uthe lifetimes of the states populated in the 
Auger decay of core ionised Mg^+(Il20)6 are given. The ICD rates are of an order of magni¬ 
tude higher than the ETMD rates. The states decaying by ETMD demonstrate comparable 
ETMD(2) and ETMD(3) efficiencies, which can be explained by the size of the system. Al¬ 
though ETMD (2) is faster in a subsystem with two neighbours, the quadratic scaling of the 
number of ETMD(3) channels with the number of neighbours n strongly enhances ETMD(3) 
for n=6. Compared to the Mg"‘+(2p“^[^D,^ S]) states, the higher lying Mg'^+(2s“^2p“^pP]) 
state exhibits a more complex ETMD pattern. Both 2s and 2p holes on magnesium may act 
as the electron acceptor, producing either Mg^'''(2p“^pP]) or Mg^’''(2s“^pS']) cations. The 
electron transfer into the 2p hole is favoured because of the greater overlap with orbitals of 
water. 

The additive approach was furthemore applied to calculate the decay rate of the core ICD- 
like processes in the Mg^+(H20)6 cluster. The core ionised Mg^+(ls“^pS'])(Il20)6 state 
decays not only locally by the Auger mechanism but is also coupled to a very large number 
of hnal states of interatomic electronic decay (e.g. core ICD-like processes). Therefore, an 
explicit computation of total and partial rates of decay for this state was not feasible. The 
rate of the core ICD-like processes was calculated for a subsystem with one water neighbour 
and extrapolated to n = 6. The Auger rate in the Mg^+(Il20)6 cluster was expected to 
be only weakly modihed compared to the atomic value, which was computed by the Fano- 
Stieltjes-ADC method. 
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III. CASCADE OF ELECTRONIC DECAY PROCESSES 


The separate electronic decay steps considered in this work can be combined into a decay 
cascade. The accumulation of the positive charge in the course of the decay cascade will have 
an impact on the electron energies and decay lifetimes. In general, the latter are expected 
to increase due to the fact that less ionisable water neighbours are available after each step 
of the cascade. The additive approach allows for a recalculation of the decay lifetime for the 
given number of ionisable neighbours n and a given decay step. Since any electronic decay 
step leads to the increase of the positive charge of the cluster by one, the hnal states of the 
decay will be destabilised relative to the initial ones. Thus, the number of available decay 
channels will further decrease and the electron energies will be shifted to lower values. 

We took this effect into account replacing ionised neighbours by point charges in the ab initio 
calculations of the energies of the final and initial states of decay. The point charge model 
is a crude approximation, however, it contains the leading term in the multipole expansion 
of the charge distribution on the ionised water molecules. Thus, the main contribution of 
the electrostatic interaction with the ionised water neighbours is included. Such calculations 
show that the last step of the most prominent pathway in the decay cascade: 

Mg^+{ls-^[^S]){H20)e Mg^+{2p-^[^D]){H20)e 

Mg^+{2p-^[^P]){H^O^UH^O), ( 12 ) 

Mg^+[^S]{H^O+UH^O)2 

is energetically possible without removing ionised water molecules produced in the first 
ETMD(3) step. 

Using the point charge model to obtain the decay widths of the individual steps (see Fig. 2 
in the main text), we estimated the time necessary for the completion of the cascade in the 
Mg3+(i_5-i)(H20)6 cluster. An Mg^+ population of 0.9 is found at 220 fs, whereas the true 
time of the completion is expected to be shorter. Fast repulsive nuclear dynamics following 
ICD and FTMD, which we did not take into account, will efficiently remove the ionised 
neighbours and open closed decay channels. 
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IV. BASIS SETS FOR ELECTRONIC STRUCTURE CALCULATIONS 


A characteristic feature of the electronic ground states of the Mg^+(H20)„ clusters as well 
as of the decaying states is the positive charge at the metal site. In this case cc-pCVTZ 
basis set on magnesium and daug-cc-pVTZ basis set on oxygen and hydrogen atoms were 
utilised for geometry optimisations and calculations of energies of the decaying states SQ. 


The calculation of the core ionisation potential of Mg^’''(H20)6 required an uncontraction 
of the atomic basis set on magnesium. In the hnal states of the electronic decay additional 
positive charge on the neighbouring water molecules is created, therefore cc-pCVTZ basis 
set for oxygen and cc-pVTZ for hydrogen atoms were chosen [ll|, Il2 |. 

Decay rate calculations required an augmentation of the standard atomic basis sets by 
Kaufmann-Baumeister-Jungen (KBJ) continuumlike functions to improve the description 
of the pseudocontinuum [1^. In particular, cc-pCVTZ (magnesium and oxygen) and cc- 
pVTZ atomic basis set (hydrogen) were augmented by Islpldlf KBJ basis functions on 
magnesium, AsApAdlf on oxygen and 2s2p2d on hydrogen. Additionally, 2s2p2d KBJ basis 
functions were located between the magnesium and oxygen atoms. An uncontraction of the 
atomic cc-pCV5Z (magnesium, oxygen) and cc-pV5Z (hydrogen) basis sets was done for the 
computation of the Auger rate, in order to describe the fast outgoing Auger electron BQ. 
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